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Abstract. Light curves of the semi-detached eclipsing bi- 
nary system TZ Eridani in the Geneva 7-colour photomet- 
ric system were analysed using the Wilson-Devinney pro- 
gramme. The physical and orbital parameters have been 
determined through a self-consistent simultaneous solu- 
tion of the seven light curves and of the radial veloc- 
ity curves of both components. The following absolute 
elements of the components are, for the primary (mass 
gainer). Mi = 1.97 ± 0.06 Mq, Ri = 1.69 ± 0.03 Rq, 
Mboh = 2.36 ± 0.09, Teffi = 7770 ± 100 K, and for the 
secondary (mass loser), M2 — 0.37 ± 0.01 Mq, R2 = 
2.60 ± 0.04 Rq, Mbob = 3.74 ±0.13, Teff, = 4570 ± 100 K. 
The semi-major axis A of the relative orbit is 10.57±0.16 
R©. The regular increase of the period is described. The 
spectral type of the components are A5/6 V (primary) 
and KO/1 III. The secondary has exhibited a long-term 
luminosity increase of about 0.06 in V between Dec 1983 
and Dec 1996. The distance to TZ Eri is evaluated to 270 
± 12 pc. 

Key words: Stars : individual : TZ Eridani - Binaries : 
eclipsing - Stars : fundamental parameters - Techniques : 
photometric - Techniques : radial velocities 



1. Introduction 

TZ Eridani (BD -6°880) is an Algol-type EA eclipsing bi- 
nary of short period, P — 2.606 days. The total primary 
eclipse is 2.5 mag deep in the V band. The estimation of 
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the spectral type of the two components is F (Cannon, 
1934), F8 (Brancewicz & Dvorak, 1980) for the primary, 
and KOIV (Kaitchuck & Park, 1988), KO (Yoon et al., 
1994) for the secondary. However, as shown in this pa- 
per, the most probable spectral types are A5/6 V for the 
primary and KO/1 III for the secondary. 

A first estimate of the physical parameters of the com- 
ponents is given by Brancewicz & Dvorak (1980), in their 
"Catalogue of Parameters for Eclipsing Binaries" , in par- 
ticular i?i — 1.51 R0 (primary component) and R2 = 1.92 
Rq. Recall that these values have not been obtained by 
a simultaneous treatment of the complete light and radial 
velocity curves. It is the aim of this paper to perform such 
a new analysis, by using the method and the computer 
programme of Wilson & Devinney (1971). 

No complete light or radial velocity curves have been 
yet published for TZ Eri. For that reason, this star was 
measured intensively in the 7-colour Geneva photometric 
system (Golay, 1980; Rufener, 1988) using the Swiss tele- 
scope at La Silla (European Southern Observatory, Chile) 
equipped with the two-channel aperture photometer P7 
(Burnet & Rufener, 1979). Moreover, the radial veloc- 
ity curve of the primary component has been determined 
with the spectrovelocimeter Coravel installed on the 
1.54 m Danish telescope at La Silla, and with the Illinois 
Cassegrain ("white") spectrograph attached to the 1 m 
Illinois reflector at Mount Laguna Observatory. In addi- 
tion, one spectrum of the secondary component has been 
obtained by using the NTT 3.5 m telescope at La Silla, 
allowing the determination of the masses of the two com- 
ponents of TZ Eri. 

TZ Eri was discovered to have a disk (Kaitchuck & 
Honeycutt, 1982) on the basis of emission in H/3 and 
lines observed during the eclipse. Kaitchuck & Park (1988) 
showed that this disk belongs to the class of the transient 
accretion disks which are produced by a collision of the gas 
stream with the mass- gaining star (primary component). 
The radius of the disk was measured by the duration of 
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the presence of H/3 line in emission after the start of the 
primary eclipse (second contact) or before the end of the 
same eclipse (third contact). The disk extension is variable 
from one eclipse to another, and in most cases between the 
trailing and leading sides. During the 12 eclipses studied, 
the disk extension r/Ri varied from 1.0 to 1.64. 

In this paper, the variation of the period of TZ Eri 
will be analysed, the variability of the components will be 
examined, and the physical parameters of the two com- 
ponents will be determined from the analysis of the light 
and radial velocity curves. 

2. Period 

The orbital period hsted in the GCVS (Kholopov, 1985) 
is P = 2.6060653 days. However, we are dealing here with 
a semi-detached system which, by definition, has a vari- 
able period due to mass transfer. During most photomet- 
ric observations done in 1984-1987 by Dr. Zdenek Kviz, 
the period remained fairly constant, since it was possible 
to obtain a well-defined lightcurve with P = 2.6061082 
days. This value was obtained using the reciprocal 9i 
test of Renson (1978) - see also Manfroid et al. (1991) 
- which had been defined above all for Ap stars (which 
vary with small amplitudes) but proved very efficient for 
a precise determination of the period of eclipsing binaries. 
The ephemeris we have adopted is : 



HJD(MinI)= (2 446109.6922 ±0.0010) (1) 
+(2.6061082 ± 0.0000020) x E 

The orbital period has varied, as shown by the times of 
minima registered over decades by amateur astronomers. 
The O—C values so obtained, and published in the BSAG 
Bulletin (e.g. Locher 1997) are shown in Figure 1 together 
with a parabolic fit which gives : 



0-C = -(0.0050 ± 0.0014) 

-(1.008 ± 0.366) X 10"'^(HJD - to) 
-F(2.433 ± 0.H3 X 10-^(HJD - t^f 



(2) 



where to = 2 446 109.6922 . The typical error on the epochs 
of primary minima in Fig. 1 is a few minutes (0.002 - 0.005 
d.). Therefore the fit only represents a mean trend, upon 
which arc superposed sudden period changes that cannot 
be accounted for by measurement errors (see especially 
the very steep rise just before JD 2450000). 

Although the period changes arc interesting by them- 
selves, they are rather a nuisance in our context, because 
the radial-velocity measurements, which were made re- 
cently, have to be put in phase with the photometric 
ones, which are much older. For this reason we have not 
used the above formula for O — C, but we have simply 
used additional photometric measurements kindly made 



0.06 



0.04 



b 0.02 



-0.02 



I I I I I , I I I , I I I I I r- 
HJD(Mm I) = 2446109.6922 + 2.6061082 E 




2000 4000 6000 8000 iq* 

HJD-2440000 

Fig. 1. 0-C diagram of TZ Eri from observations made by 
amateurs and published in the BSAG Bulletin (two points are 
defined by Geneva photometry). A parabolic fit is superposed, 
showing the regular increase in the period. Short arrows define 
the intervals of photometric measurements; long arrows define 
the interval of radial-velocity observations. 



by Marc Kiinzli in November 1996 with the same equip- 
ment. He has made 36 new multicolour measurements, 
several of which during the primary minimum. Using the 
code EB0P16 (Etzel, 1989) and the adopted ephemeris, 
we adjusted the A9 parameter (phase correction for the 
position of the primary minimum) for these recent data, 
as well as for the old data alone. The difference is : 



A(j) 



-0.01734 ±0.00020 



(3) 



and represents the phase correction to be applied to the 
1996 data, to fit them into the adopted ephemeris. Al- 
though the Vr observations have been made one year ear- 
lier than the new photometric ones, we neglect the slight 
period change that may have occurred in between, com- 
pared to the change that has taken place between the old 
photometric measurements and the Vr ones. Therefore, 
the above phase shift was applied as such to the Vr data. 

3. Photometric data and variability of the compo- 
nents 

Geneva 7-colour photometric measurements of TZ Eridani 
were obtained from Dec 10, 1983 to Dec 10, 1996, using 
the Swiss 70 cm telescope at the European Southern Ob- 
servatory (ESO), La Silla, Chile. During this period, 393 
measurements of weight ^ > 1 have been obtained (see 
Rufener, 1988, for the definition of the weight q). These 
data are listed in Table 1, together with the 36 additional 
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photometric measurements obtained in Nov-Dec 1996 (see 
Section 2). 

The magnitudes in each of the seven filters are ob- 
tained from the visual magnitude V and the six colour 
indices in the following manner : 

i = V -[V - B] + [i- B] (4) 

with i representing one of the seven filters U, B, V, Bi, 
B2,Vi, G. Remember that the Geneva [U-B] and [B-V] 
indices are not normalized to zero for an AOV star as it 
is the case for the Johnson UBV indices. It is possible to 
calculate the magnitude of the primary (mass gainer) by 
subtracting the flux of the secondary (mass loser), at the 
bottom of the primary eclipse, from the flux of the both 
components measured together (outside the eclipses). This 
calculation has been made for each of the seven Geneva 
magnitudes. In order to minimize the effects of a possible 
long-term variability (see the end of this Section), only the 
data obtained in December 1983 and January 1984 have 
been used. The results are given in Table 2. 

It is interesting to compare the observed uncertain- 
ties with the mean precision of the measurements made 
in Geneva photometry. Rufcncr (1988, Fig. 2) has shown 
the shape of the mean relation cry vs. V obtained for the 
non- variable stars, in particular the progressive increase 
of ay with increasing V, for stars fainter than V ^ 9. The 
same relation can be applied to the seven Geneva magni- 
tudes. Fig. 2 of this paper shows a new calculation of this 
relation, based on the up to date version of our photomet- 
ric database. On the same figure are plotted the observed 
values for TZ Eri given in Table 2. A correction has been 
applied to the uncertainties of the magnitudes for the sec- 
ondary, because the measurements obtained during the to- 
tality of the primary eclipse had shorter integration time 
than the other ones (4 minutes instead of 12 minutes). The 
conclusion is that the uncertainties on the measurements 
of TZ Eri are in agreement with the expected precision. 
Thus, the components do not exhibit a short-term vari- 
ability, i.e on a time-scale shorter than about 20 orbital 
periods. 



Table 2. The seven Geneva apparent magnitudes of TZ Eri 
(measures of December 1983 and January 1984). 



Mag. 


TZ Eri A+B 


Secondary 


Primary 




Observed 


Observed 


Calculated 


U 


10.540 ± 0.013 


14.87 ±0.17 


10.560 ± 0.016 


Bl 


9.987 ±0.009 


13.980 ±0.050 


10.015 ±0.011 


B 


9.011 ± 0.009 


12.584 ± 0.031 


9.052 ±0.010 


B2 


10.409 ± 0.007 


13.670 ± 0.035 


10.464 ± 0.009 


VI 


10.299 ±0.011 


12.809 ±0.021 


10.406 ±0.014 


V 


9.576 ±0.010 


12.067 ±0.019 


9.691 ±0.013 


G 


10.679 ± 0.016 


12.960 ± 0.018 


10.821 ± 0.021 




_l I I I I I I I I I I I L_ 

10 12 14 16 

mag. 



Fig. 2. Variation of the mean precision a with the magnitude, 
in the case of Geneva photometric measurements. The solid 
thick line refers to the mean value and the solid thin lines to 
the 1 s.d. level. Dots and small squares concern respectively 
the seven magnitudes of TZ Eri A±B outside eclipses and of 
the secondary (during the totality of the primary eclipse) . Big 
squares represent the estimated values of a for the secondary 
which ought to have resulted from a "normal' integration time 
of the measurements, i.e. 12 minutes (normal) instead of 4 
minutes (during the primary eclipse) . 

The long-term photometric behaviour of both compo- 
nents has been analysed by comparing our photometric 
data obtained at 4 epochs, corresponding to the inten- 
sive monitoring of the eclipses : Dec 1983 to Jan 1984, 
Jan 1985, Nov 1987 and Nov-Dec 1996. Table 3 gives the 
mean values of V, [B — V] and \U — B] at each of these 
epochs for the both components. It appears that : 

1. The secondary exhibited a long-term luminosity in- 
crease (0.06 in V) between Dec 1983 and Dec 1996. 
The colour variations are large, especially in \U — B], 
but not significant due to the large standard deviation. 

2. The primary did not show any long-term variation in 
magnitude or in colours. 

The observed variations of TZ Eri secondary are sim- 
ilar to those studied by Olson & Etzel (1993) in six cool 
subgiant secondaries of totally eclipsing Algol systems. 
They noted that the fluctuations increase with decreasing 
orbital period, or with increasing rotational velocity, sug- 
gesting that rotationally induced magnetic activity could 
be the origin of these brightness variations. 

In order to minimize the effects of the long-term vari- 
ation of the secondary luminosity, only the data obtained 
before HJD 2447200 (Feb 1988) have been used for the 
eclipse analysis (see Section 6). 
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Table 3. Long-term behaviour of the components of TZ Eri in V magnitude and Geneva colours [B — V] and [U — B]. Only 
the variation of the secondary in V is significant. 





Secondary (mass loser) 
V [B-V] [U- B] 


Primary (mass gainer) 
V [B-V] [U- B] 


Dec 1983 

Jan 1985 
Nov 1987 
Dec 1996 


12.067 0.523 2.281 
±0.019 ±0.046 ±0.175 
12.045 0.466 2.184 
12.035 0.453 1.990 
12.007 0.477 2.189 


9.691 -0.639 1.508 
±0.013 ±0.023 ±0.026 

9.691 -0.639 1.561 

9.692 -0.637 1.561 
9.695 -0.635 1.547 



4. Classification of the components and interstellcir 
extinction 

A photometric classification of the components of TZ Eri 
can be obtained because the primary eclipse is total, by 
using various calibrations based on the colour indices and 
reddening free parameters of the Geneva photometric sys- 
tem. 

From the magnitudes in Table 2, the various Geneva 
colour indices and reddening free parameters d, A and g 
(see Golay, 1980) have been calculated. Then, the tech- 
nique of the photometric boxes (Golay et al., 1969; Nico- 
let, 1981a) has been applied to determine the intrinsic 
characteristics of the primary component. Recall that this 
technique is specific to the Geneva photometric system 
and depends on the homogeneity and precision of its mea- 
surements. The assumption in this technique is that the 
properties of the stars lying within a specified photomet- 
ric neighbourhood (see Nicolet, 1994) can be equated with 
each other, provided that the radius of the photometric 
box is small enough. We applied the technique to the pa- 
rameters d, A and g of the primary and found in the entire 
photometric database 40 twin stars, i.e. stars having very 
similar values of the 3 parameters, the radius of the pho- 
tometric box having been chosen equal to 0.015 mag. It 
is assumed that these stars are intrinsically similar to the 
primary of TZ Eri. The same method was applied to the 
classification of the components of the eclipsing RS CVn- 
type system RZ Eridani (Burki et al., 1992). 

It is especially noteworthy that 6 of these twin stars 
belong to open clusters : aPcr, Hyadcs, Pracscpc, Pleiades 
(2 stars) and NGC 6405. This allowed to determine unan- 
biguously the intrinsic colours of these stars, since these 
clusters have well determined distances and interstellar 
extinctions. By using the values given by Nicolet (1981b), 
we derived the intrinsic colours of the twin stars, and thus 
also of the primary of TZ Eri, in particular [B2 — V1]q — 
0.011 ± 0.016. According to the relations between spec- 
tral types and intrinsic Geneva colours or parameters by 
Hauck (1994) and to the calibration of Geneva photom- 
etry based on Kurucz's atmosphere models by Kiinzli 
et al. (1997), the estimated spectral type of the pri- 
mary is A5/6 V, its effective temperature is 7770 ± 100 K, 



its gravity is log 5 = 4.40 ± 0.07 and its metallicity is 
[M/H] = 0.13 ± 0.09. On the other hand, the colour ex- 
cess is E[B2 - VI] = 0.047 ± 0.017. 

For the secondary, we have calculated the intrinsic 
colours by using the measured values (see Table 1) and 
the colour excesses obtained for the primary. We obtained 
in particular [B2 — Vl]o = 0.756, an estimated spectral 
type of KO/1 III and an effective temperature of 4520 K. 

5. Radial velocity curves 

The primary's radial velocity curve has been measured 

with the CORAVEL scanner (Baranne et al. 1979) in 
December- January 1995-1996. The cross-correlation dip is 
rather wide and shallow due to the fast axial rotation, so 
that the scatter of the residuals around the fitted curve is 
rather large (several kms^^). The secondary component 
remained invisible to CoRAVEL, due to its low luminos- 
ity in the blue. About one year later, 14 new observations 
have been done by E.G. Olson using the Illinois Casscgrain 
("white") spectrograph : 38 cm focal length camera, 831 
line mm^^ grating, Texas Instruments 800x800 CCD de- 
tector cooled with liquid nitrogen; the CCD was operated 
by a Photometries Ltd. controller; spectral resolution in 
the first order is 0.04Aper pixel; Reductions were done 
with IRAF. 

There has been a problem in fixing the zero-point of 
the radial velocities obtained with this instrument, so that 
we fitted an orbit separately to the corresponding data to 
obtain the apparent systemic velocity. Then we applied to 
these Vr values a uniform shift equal to the difference be- 
tween this apparent systemic velocity and that obtained 
with CoRAVEL measurements alone. Indeed, we are con- 
fident in the Coravel Vr scale, because several standard 
stars have been observed each night and the small in- 
strumental drifts (±1 kms~^ at most) are well controlled. 
Therefore, the uncertainties quoted in Table 5 may be 
slightly optimistic, because they refer to an orbital solu- 
tion which assumed a perfect correction to Olson's radial 
velocities. 

In order to see the secondary star's spectrum and to 
obtain the mass ratio of the components, we asked Dr. Di- 
dier Raboud to observe TZ Eri in the vicinity of a quadra- 
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ture with the NTT telescope at ESO. He could indeed 

take one spectrum, with an exposure time of 10 min, on 
11th November 1995, using the EMMI spectrograph in the 
REMD mode, with Grism #5, Grating #10 and a slit mea- 
suring 1" X 6" ; in this configuration, the resolving power 
is i? = 28000 and the wavelengths range between 4013 
and 6606 A. The detector was CCD #36 (ESO numerota- 
tion), a thin, back-illuminated Tektronix TK2048EB chip 
with 24 X 24 /xm pixels. The spectrum has been reduced 
at Geneva Observatory by Mr. Michel Studer, using the 
Tacos software developed by Dr. Didier Queloz for the 
Elodie spectrograph at Observatoire de Haute-Provence. 
The radial velocities were obtained by cross-correlation 
between the observed spectrum and a binary mask opti- 
mized for FO-type stars, which yielded two dips, one for 
each component. Thanks to the long wavelength interval 
extending well into the red, the cool companion is easily 
seen in the correlation function. The dips are only 3.6% 
and 1.1% deep for the primary and secondary respectively, 
but the S /N ratio of the correlation fimction is better than 
500. A KO-type mask was tried as well and yielded the 
same result, but with no improvement. 

The journal of the radial velocity observations is given 
in Table 4 and the radial velocity curve is shown in Fig- 
ure 3. Additional spectra would of course be welcome to 
complete the secondary's Vr curve, but the single point we 
have suffices to constrain the mass ratio to a precision of 
about 7% (g = 0.193 ± 0.013). The orbital elements are 
given in Table 5. 




0.5 1 

phase 

Fig. 3. Radial velocity curve of TZ Eridani. The black dots 
(CORAVEL observations) and the plus signs (Olson's observa- 
tions) represent the primary component, while the open dot 
represents the secondary. The phases are those of Table 4. 



Table 4. Journal of the radial velocity observations of TZ Eri. 
The phases are computed from the ephemeris given by Equa- 
tion 1 and corrected for period change according to Equation 3. 



rlju 


comp. 


Vr±a 


Instrument 


phase 


— z 4UU UUU 




[km s 




or observer 




50033.785 


A 


18.08 


2. 


Emmi 


0.712 




B 


-173.40 


3. 






50039.606 


A 


0.67 


2.87 


CORAVEL 


0.945 


50042.599 


A 


-27.19 


3.08 




0.094 


50083.578 


A 


10.32 


2.50 




0.818 


50084.542 


A 


-42.17 


2.55 




0.188 


50084.708 


A 


-48.40 


2.89 




0.252 


50085.544 


A 


-1.72 


1.87 




0.572 


50085.707 


A 


11.17 


2.30 




0.635 


50086.550 


A 


3.38 


1.56 




0.958 


50087.540 


A 


-41.39 


2.17 




0.338 


50087.684 


A 


-38.53 


2.46 




0.394 


50088.551 


A 


17.60 


1.96 




0.726 


50088.709 


A 


15.04 


2.63 




0.787 


50097.630 


A 


-42.39 


2.23 




0.210 


50098.606 


A 


3.24 


2.37 




0.584 


50110.324 


A 


-29.78 


3.29 




0.081 


50407.863 


A 


-47.05 


2.5 


E.G. Olson 


0.251 


50407.880 


A 


-44.75 


2.5 




0.257 


50407.897 


A 


-46.65 


2.5 




0.264 


50408.730 


A 


^ 4^ 


9 n 




0.583 


50411.717 


A 


18.55 


2.0 




0.729 


50411.745 


A 


20.75 


2.0 




0.740 


50411.872 


A 


16.35 


2.0 




0.789 


50412.697 


A 


-34.05 


2.0 




0.106 


50412.719 


A 


-35.95 


2.0 




0.114 


50412.853 


A 


-40.95 


2.0 




0.165 


50413.735 


A 


-15.05 


2.5 




0.504 


50413.931 


A 


3.15 


1.5 




0.579 


50415.700 


A 


-39.65 


2.5 




0.258 


50415.725 


A 


-40.15 


2.5 




0.267 



6. Rotation and Metallicity 

The rotational velocity of the primary deduced from the 
standard calibration of the width of the correlation dip 
measured with Coravel (Benz & Mayor, 1984), is 

t;sini = 38.1±2.2A;ms"^ (5) 

which is slightly faster than the synchronous rotation 
value, since with a radius R = 1.7Oi?0 (see Section 6), the 
expected equatorial velocity of the primary is 33.0 kms~^. 
It is not exceptional for Algol gainers to rotate faster than 
synchronously as shown by Koch et al. (1965). However, 
the difference being only 2.3a, we cannot firmly conclude 
that a departure from synchronism exists in TZ Eri. On 
the other hand, v sin i may be overestimated, due to sev- 
eral factors. First, the rotational velocity is so large that 
we hardly reach the continuum of the correlation func- 



6 



F. Barblan et al.: Physical parameters of the Algol system TZ Eridani 



Table 5. Orbital elements of the binary. For each component, the second line gives the estimated standard deviations of the 
parameters. A null uncertainty means that the corresponding parameter has been fixed before the convergence. The period has 
been fixed to the value it had in the middle of the time spanned by the Vr observations. 



Star name 


P 


To [HJD 


e 


Vo 






Ml, 2 sin'^i 


ai,2 sin i 


N 


(0-C) 




[days] 


-2450 000] 




[km s-i] 


n 


[km s-i] 


[Mo] 


[10'^ km] 




[km s"^] 


TZ Eri A 


2.606132 


83.4009 


0.000 


-13.10 


0.0 


31.93 


1.741 


1.144 


30 


3.04 




0.0000000 


0.0000 


0.000 


0.55 


0.0 


0.70 


0.087 


0.025 






TZ Eri B 












165.32 


0.336 


5.92 


1 
















3.73 


0.017 


0.13 







Table 6. Characteristics of the 7 Geneva passbands (Rufener 
& Nicolet, 1988). Ao is the mean wavelength and fi is the sec- 
ond-order moment (approxiniatoly half the passband width). 



Passband 


Ao[A] 


M[A] 


U 


3464 


159 


Bl 


4015 


188 


B 


4227 


282 


B2 


4476 


163 


VI 


5395 


202 


V 


5488 


296 


G 


5807 


200 



Table 7. Adopted values of the logarithmic limb-darkening 
parameters which were kept fixed in the least-squares solution. 

For the primary, only the y parameters were fixed (the fitted 
X parameters are listed in Table 8). For the secondary, both x 
and y parameters had to be fixed. 



Passband 


2/1 


X2 


3/2 


U 


0.243 


0.940 


-.515 


Bl 


0.303 


0.863 


-.240 


B 


0.303 


0.858 


-.240 


B2 


0.303 


0.868 


-.240 


VI 


0.263 


0.847 


-.043 


V 


0.263 


0.825 


-.043 


G 


0.263 


0.819 


-.043 



tion, which can bias the fit. Second, the profile is fitted 
by a gaussian, while the rotational profile itself does not 
have this shape, even if it is convoluted by a gaussian in- 
strumental profile (the calibration was intended for slow 
rotators, with vsini generally slower than ~ 20 kms~^). 
Third, the calibration has been devised above all for cool, 
solar-type stars (although it does include a temperature 
term) and we apply it to a late A star, i.e. at the very 
border of its validity. Fourth, the vsini determination is 



based on the hypothesis of a solar-type macroturbulence; 
this is probably not valid in a semi-detached short period 
system and, thus the value given in (5) could be an over- 
estimate. 

We did not attempt to derive the projected rotational 
velocity of the secondary component by using the corre- 
lation dip obtained with the NTT spectrum, because a 
new calibration would have been needed which is beyond 
the scope of this paper, and the correlation dip of the 
secondary is in any case very shallow. 

It is interesting to notice that the "surface" or equiv- 
alent width of the CORAVEL autocorrelation dip of the 
primary is the same as for single late-A stars. We have 
W = 2.31 ± 0.14 kms^\ while HD 2628, type A7III and 
[Fe/H]=-0.02, has W = 2.41 and HD 110379, type FOV 
and [Fe/H]=-0.07, has W = 2.49. This implies that the 
mctallicity of TZ Eri is close to the solar value, in agree- 
ment with the result obtained in Section 4 on the basis of 
the photometric analysis. 

7. Photometric solution 

The photometric solution for TZ Eri eclipses was obtained 
with the Wilson-Dcvinncy WD programme (Wilson & 
Dcvinncy, 1971; Wilson, 1992), using the version revised 
in 1995. We used the WD programme in Mode 5, adapted 
for semi-detached systems , allowing a simultaneous com- 
putation on the seven Geneva photometric light curves, 
based on 393 measurements (see Section 3), and the ra- 
dial velocity curves of both components (see Table 4). 

Lobe filling of the cool loser component was assumed. 
The semi-major axis of the relative orbit was initially set 
to ^ = 10.2 Rq as calculated from radial velocity re- 
sults (see Section 5). Orbital eccentricity is fixed to zero 
and longitude of the periastron to 90° . Primary star tem- 
perature was fixed to 7770 K as determined in Section 4. 
For both components, the stellar atmosphere models of 
Kurucz (1994) integrated by Nicolet (1998) through the 
Geneva photometry passbands (Rufener & Nicolet, 1988) 
have been used. 
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Table 8. Adjusted parameters of the system TZ Eri from the Wilson-Devinney programme. The parameters are : 
the temperatures Ti and T2 of the primary (hot) and secondary (cool) components, 
the semi-major axis A of the relative orbit, 
the orbital inclination i, 

the mass ratio q (the uncertainty assumes exact lobe filhng for the secondary), 

the potential of the surface of the primary component f2i (for the units, see Section 7 and WD programme), 
the potential of the surface of the secondary component Q.2 (for the units, see Section 7 and WD programme), 
the exponent of the gravity darkening law 52, 

the normalized monochromatic luminosity in the seven Geneva passbands for the primary Li/(Li + L2), 
the normalized monochromatic luminosity in the seven Geneva passbands for the secondary L^jiL^ + L2), 
the center-to-limb darkening factors for the primary x\ in the seven passbands. 



Parameters 


Values 

(Ti = 7770is:) 


Uncertainties 


Values Differences 
(Ti = 7670Ji:) "7670-7770" 


Values 

(Ti = 7870Ji:) 


Differences 
"7870-7770" 


Ti [K] 


7770 




7670 


-100 


7870 


+100 


J 2 

A [i?0] 

q = M2/M1 


40Do 

10.57 

86.73 
0.1865 


4-0 

±0.16 
±0.03 
±0.0003 


4614 
10.59 
86.77 
0.1862 


51 
0.02 
0.04 
-0.0003 


4DDy 

10.56 
86.71 
0.1867 


lUD 

-0.01 
-0.02 
0.0002 


Qi 
Q2 
92 


6.46 
2.20 
0.322 


±0.01 


6.45 
2.20 


-0.01 
0.0 


6.49 
2.20 


0.03 
0.0 


Li/{Li+L2) 
U 
Bl 
B 
B2 
VI 
V 
G 


0.9687 
0.9644 
0.9481 
0.9312 
0.8763 
0.8676 
0.8420 


±0.0010 
±0.0003 
±0.0003 
±0.0003 
±0.0004 
±0.0005 
±0.0003 


0.9692 
0.9644 
0.9481 
0.9313 
0.8765 
0.8678 
0.8422 


0.0005 

0.0 

0.0 

0.0001 
0.0002 

0.0002 
0.0002 


0.9692 
0.9644 
0.9481 
0.9313 
0.8765 
0.8679 
0.8421 


0.0005 

0.0 

0.0 

0.0001 
0.0002 

0.0003 
0.0001 


L2/{Li+L2) 
U 
Bl 
B 
B2 
VI 
V 
G 


0.0313 
0.0356 
0.0519 
0.0688 
0.1237 
0.1324 
0.1580 


±0.0023 
±0.0006 
±0.0005 
±0.0005 
±0.0006 
±0.0006 
±0.0004 


0.0308 
0.0356 
0.0519 

0.0687 
0.1235 
0.1322 

0.1578 


-0.0005 

0.0 

0.0 

-0.0001 
-0.0002 
-0.0002 

-0.0002 


0.0308 
0.0356 
0.0519 

0.0687 
0.1235 
0.1321 

0.1579 


-0.0005 

0.0 

0.0 

-0.0001 
-0.0002 
-0.0003 

-0.0001 


XI 

U 

Bl 

B 

B2 

VI 

V 

G 


0.655 
0.770 
0.745 
0.735 
0.533 
0.542 
0.463 


±0.037 
±0.009 
±0.009 
±0.009 
±0.013 
±0.015 
±0.012 


Determinant of normal equations: 1.23 • 10 ^® 

Predicted sum of weighted residuals : 0.00639 
Observed sum of weighted residuals : 0.00640 



The bolomctric albedos for hot and cool components 
were taken at the theoretical value of 1.0 and 0.5 respec- 
tively (radiative and convective cases). Noise was set to 1 
because scintillation can been neglected, the photon noise 
being dominant (see Bartholdi et al., 1984). The grid reso- 



lution values were taken as 20, 20, 20, 20 for Nl, N2, NIL 
and N2L respectively (see WD programme). For each of 
the seven Geneva magnitudes, the passband mean wave- 
length Ao was the one calculated by Rufencr & Nicolet 
(1988) as shown in Table 6. We assigned the Aq value for 
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the B filter to the passband of the spectrographs used 
for the radial velocity measurements. Stellar rotation is 
assumed to be synchronized for both components. For 
both the primary and secondary components, a logarith- 
mic limb-darkening law of the form : 

/ = /o(l — a; -|- a; cos 9 — y cos 6 ln(cos 6)) (6) 

was assumed (Van Hamme, 1993). For the secondary, both 
X and y parameters were fixed to their theoretical values, 
interpolated from Table 2 of Van Hamme (1993). Indeed, 
the secondary minimum is too shallow to allow the de- 
termination of X2- For the primary, the yi parameter was 
fixed to its theoretical value, because it is not possible 
in practice to determine both xi and j/i parameters since 
they are strongly correlated. The xi parameter, on the 
contrary, was left free and was fitted in the least-squares 
procedure, for each of the seven passbands. The adopted 
values of yi , X2 and y2 are listed in Table 7, while the fitted 
values of xi are given in Table 8. The gravity darkening 
exponent 172 is relevant to the secondary component and 
is defined by the expression (Wilson & Biermann, 1976) : 

where Tiocai is the local effective temperature, T^.p. is the 
temperature at a reference point which here is the com- 
ponent's pole, and aiocab o,r.p. are the accelerations due to 
gravity at the same respective points. 

The adjustable parameters are then semi-major axis 
A, inclination i, mass ratio q = M2/M1, cool star temper- 
ature T2, hot and cool star luminosities Li and L2 in each 
passband, cool star gravity darkening exponent 92 (see 
Wilson & Biermann, 1976), hot star limb darkening coef- 
ficient xi in each passband and potential fl\ at the hot star 
surface. In reality, 17 is a non-dimentional parameter which 
is a linear function of the true potential \l/ (Kopal, 1959; 
Wilson & Devinney, 1971). First, preliminary light curves 
were generated using the Light Curve sub-programme of 
WD until a satisfactory fit to the observed light curves 
was obtained. Second, the method recommended by Van 
Hamme & Wilson (1986) was applied for the determina- 
tion of the temperature and luminosity of the cool compo- 
nent. In a first step, T2 was determined with IPB option 
(see the WD programme) equal to 0, i.e. the luminosity is 
coupled to the temperature. In a second step, temperature 
and luminosity were de-coupled (IPB = 1), and it was 
possible to determine L2 ■ Finally, the following sets of pa- 
rameters were adjusted : first, successively (T2, Qi, i), (Li, 
L2), {A, i), {xi, 52)5 and then simultaneously (A, i, Cli, 
q, Li, L2), as recommended by the users of the WD pro- 
gramme. The solution was carried out to the point where 
the probable errors of all these parameters were smaller 
than the computed parameter corrections. 

Since the photometric temperature of the primary star 
is determined with an uncertainty of ±100 K (see Sec- 
tion 4), it is interesting to calculate the solution also for 



the two Ti values 7670 K and 7870 K. The results are pre- 
sented in Table 8. The differences are very small, except 
of course on T2. 

It is especially noteworthy that, for Ti = 7770 K, the 
value of T2 (4563 ± 2 K) is very close to the value obtained 
in Section 4 on the basis of the photometric calibrations. 
The mass ratio q = M2/M1 = 0.1865 ± 0.0003 in Table 8 
is a little smaller than the value derived from the radial 
velocity analysis {K1/K2 = 0.193 ± 0.013, see Table 5), 
but still well within the uncertainty of the spectroscopic 
determination. A similar diiference can be noted on the 
semi-major axis of the relative orbit A = ai+a2 (10.16 ± 
0.22 Rq in Table 5 and 10.57±0.16 R© in Table 8). These 
differences are not surprising. We have to recall that the 
radial velocity curve for the secondary component is based 
on only one measurement. The constraints imposed by 
the photometric measurements of the eclipses produced 
an improvement of the orbital and physical parameters of 
TZ Eri. 

The small difference between the photometric and 

spectroscopic values of q could indicate that the lobe fill- 
ing by the secondary is not complete. If this were true, the 
real value of the uncertainty on q would be larger than the 
photometric value (0.0003), may be as large as the spec- 
troscopic value (0.013). Nevertheless, test solutions have 
been obtained with the cool star slightly detached from 
its Roche lobe, and the results are less satisfactory than 
those obtained in Mode 5 of the WD code : as the cool 
star became farther detached from the lobe, solution er- 
rors grew. This confirms the semi-detached status of the 
system. 

The uncertainties on the derived parameters are of 
three types : 

- Intrinsic, i.e. resulting from the mathematical anal- 
ysis of the light and radial velocity curves (e.g. ± 2 K on 
T2, see Table 8). Note that the intrinsic errors on i?i_2 for 
A fixed are an order of magnitude smaller than the values 
in Table 9. 

- Strongly correlated with the determination of A. For 
example, in Table 9, the errors on Mi_2, R1.2 and log 511^2 
almost entirely come from the uncertainty on A, and are 
pair wise strongly correlated (as well as with A). 

- Depending on the photometric determination of Tes^ . 
For example, the values of Mboii,2 strongly depend on the 

adopted value for Toff ^ . 

The light curves for the seven filters U, Bi, B, B2, 
Vi, V and G are shown in Figs. 4 and 5. The quality of 
the fits is clearly extremely good. Figs. 6 and 7 present 
two views of TZ Eri. Fig. 6 is a "classical" representation 
of the two components in the equatorial plane. The 3- 
dimentional representation of the potential (Fig. 7) allows 
a better understanding of the meaning of the Lagrange 
points and of the path for the flow of the material from 
the secondary when it fills its Roche lobe. 
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Table 9. Computed parameters of the system TZ Eri from the 
Wilson-Devinney programme. The parameters are, for the hot 
primary (1) and the cool secondary (2) components : the mass 
M, the mean radius R, the surface gravity log the bolometric 
magnitude Mboi and the various radii (in units of semi-major 
axis) of the deformed components, i.e. pole (perpendicular to 
the orbital plane), point (in the direction of the other compo- 
nent), side (in the orbital plane, in the direction perpendicular 
to the direction of the other component) and hack (in the di- 
rection opposite to the other component). The uncertainty on 
Mboi is calculated with an uncertainty on of ± 100 K. 



Parameters Values Values Values 

(Ti = (Ti = (Ti = 

7770^:) 7Q70K) 787(iK) 



Ml [M 


0] 


1 


97 ± 


0.06 


1 


98 


1 


97 









37 ± 


0.01 





37 





37 


Ri [Rq] 


1 


69 ± 


0.03 


1 


69 


1 


68 


R2 




2 


60 ± 


0.04 


2 


60 


2 


60 


log 51 




4 


28 ± 


0.03 


4 


28 


4 


28 


log 52 




3 


17 ± 


0.03 


3 


17 


3 


17 


Mboii 




2 


36 ± 


0.09 


2 


41 


2 


32 


Mbol2 




3 


74 ± 


0.13 


3 


69 


3 


64 


'^polei 


[A] 





1592 


± 0.0003 





1595 





1585 


'^pointi 







1598 


± 0.0003 





1601 





1591 


^sidei 







1596 


± 0.0003 





1599 





1589 


^backi 







1598 


± 0.0003 





1601 





1591 


^^pole2 


[A] 





2282 


± 0.0001 





2281 





2283 


^^point2 







34 


± 0.01 





3349 





3351 


^8ide2 







2374 


± 0.0001 





2373 





2374 


^back2 







2696 


± 0.0001 





2695 





2697 



8. Discussion 

8.1. Limb- darkening 

Our adjusted values of xi (sec Table 8) compare relatively 
well with the theoretical ones, in the sense that the dif- 
ference hardly exceeds 0.1, or 15-20% of the value (except 
for V where the difference reaches 0.15). Nevertheless, our 
values are systematically smaller than the theoretical ones 
(computed for a solar chemical composition). The discrep- 
ancies are larger than the formal errors indicated in Ta- 
ble 8, which may partly, but probably not entirely be ex- 
plained by the fact that the errors listed should represent 
only lower limits to the true ones. 

Fig. 8 illustrates the dependance of xi on wavelength 
for a theoretical star with T^s = 7770 K and \ogg = 4.28, 
for the passbands of Stromgren's uvby system and John- 
son's UBV system (the values have been linearly inter- 
polated in Table 2 of Van Hamme, 1993). The fitted val- 



ues for the Geneva U,B1,B,B2,V1,V,G passbands are 
shown for comparison. 
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Fig. 6. Schematic view of the system in the equatorial plane 
(coordinates in Rq). Some equipotential lines and the Lar- 
grange points are indicated. The secondary star fills its Roche 
lobe. 



The theoretical dependance of xi on wavelength sat- 
isfies very well the observations from the U up to the B2 
passband, except for a small, systematic vertical shift. For 
the Vl to G passbands, however, the observed xi param- 
eter is much smaller than the theoretical one. Since the 
error bars are not larger than the symbols in Fig. 8 (ex- 
cept marginally for U), there is undoubtedly a significant 
discrepancy which remains to be explained. We obtained 
a simultaneous solution in the seven passbands, thus at 
least the relative values of xi should be reliable. But, even 
considering only differential values (e.g. xi{V) — xi{B)) a 
difference with the theoretical predictions does remain. 
Unfortunately, we did not see any recent discussion in the 
literature about the reliability of empirical limb-darkening 
coefficients in Algol-type systems. 



8.2. Distance 

From the values of V (Table 2), Mboi (Table 9), E[B2 - 
VI] and [B2 — V1]q (Section 4), we derive a distance of 
270 ± 12 pc for TZ Eri, by adopting a zero value for the 
bolometric correction BC of the primary, according to the 
BC-colour relation of Flower (1977). 
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Fig. 7. Schematic view of the gravity potential f2 (see Sec- 
tion 7 for the definition) of the TZ Eri system. The horizontal 
coordinates are in Rq. 



1 



+ 

0.6 - 

+ + 

1 

0.4 - 



0.2 - 



' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 1 

3000 4000 5000 6000 

\ [A] 

Fig. 8. Limb-darkening xi coefficient of the primary as a func- 
tion of wavelength. Points are theoretical heterochromatic val- 
ues for the uvby (open dots) and UBV passbands (full dots), 
while crosses are for empirical values determined in this work 

for the passbands of the Geneva system. Notice that the theo- 
retical points for the v and the V bands are superimposed. 

8.3. Algol-type stars 

Sarma et al. (1996) gave a detailed discussion on the evo- 
lutionary status of Algol components, on the basis of a 
comparison of their global parameters (masses, radii, lu- 
minosities, temperatures) with those of normal stars. As 
these authors noted : 

1. While primaries (mass gainers) of some semi-detached 
systems are overluminous, oversized and hotter for 



their masses, some others are under luminous, smaller 
and cooler. 

2. Except for a few cases, the primaries are lying either 
on the main sequence or near it. 

3. The secondaries (mass losers) have evolved off the 
main sequence (they were originally the more massive 
components). 

Figs. 9 and 10 present the Mass Luminosity diagram 
(logi vs. logM) and the HR diagram (logi vs. logT) for 
the Algol systems in the lists by Sarma et al. (1996) and 
Maxted & Hilditch (1996) having a mass of the primary 
smaller than 4 Mq. The 3 systems with primaries off the 
main sequence (RZ Cnc, AR Mon, AW Peg) are not repre- 
sented. The underluminous (black dots) and overluminous 
(black squares) have been separated in Fig. 10 by using 
the ZAMS (Schaller et al., 1992) to make the separation. 
We can see that : 

1. TZ Eri belongs to the group of the underluminous pri- 
maries. 

2. The mass of TZ Eri's primary is the smallest of the 
group of these underluminous primaries. 

3. The luminosity of TZ Eri's primary is that of a main 
sequence star of 1.70 Mq located near the ZAMS (see 
Fig. 10), while the determined mass is 1.98 Mq (see 
Table 9). 

4. With respect to a main sequence star of the same mass, 
TZ Eri's primary is underluminous by ~0.3 in logi (or 
fainter by ~0.75 in Mbd) and cooler by ~900 K in Tea. 

5. The mass transfer from the secondary to the primary 
has been important, taking into account the present 
masses of the components, i.e. 0.37 Mq and 1.98 Mq. 
The minimum value of this transfered mass is given 
by: 

AM = Ml - (Ml + M2)/2 ~ 0.80Mq. 
9. Conclusion 

The simultaneous adjustment of the light curve model on 
high precision photometric data in seven passbands, rang- 
ing from 3400 to 6000 A, put strong constraints on the 
physical and orbital parameters of TZ Eri. This is our first 
analysis of an Algol system on the basis of measurements 
made in the Geneva photometric system. The essential re- 
sults are presented in Figs. 4 and 5, Tables 8 and 9, and 
in the Abstract. The quality of the fits on the light curves 
and the uncertainties on the parameters show that this 
analysis has been successful. 

A further improvement on the determination of the 
physical parameters of the components of TZ Eri, in par- 
ticular the mass ratio, would result from the acquisition 
of additional high S/N spectra of the secondary cool com- 
ponent. Indeed, the present study is based on only one 
spectrum of the secondary and this is not completely sat- 
isfying. However, this spectrum was obtained using the 
3.5 m NTT telescope at La Silla and we were extremely 
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Fig. 9. Theoretical HR diagram for the components of Algol 
systems with very precisely determined physical parameters 
according to the lists by Sarma et al. (1996) and Maxtcd & 
Hilditch (1996), and to this paper for TZ Eri (see Section 8.3). 
Primaries are identified with filled symbols and secondaries 
with open symbols. TZ Eri is identified with plus symbols. 
Dots and squares refer to the systems having respectively un- 
der- and overluminous primaries with respect to their masses 
(see Fig. 10). The Zero- Age Main Sequence and 3 evolutionary 
tracks (3.0, 2.0 and 1.5 Mq) are drawn, according to Schaller 
et al (1992). 

lucky to have access to this large telescope for this single 
measurement. 

Accordind to Maxted & Hilditch (1996), the number 
of Algol-type binary systems for which the absolute pa- 
rameters are determined on the basis of a self-consistent 
solution of both the light and radial velocity curves is very 
limited. They mentioned only 9 systems for which masses, 
radii and luminosities are known to accuracies typically 
better than 5%. In this context, it is important to enlarge 
the sample of very well known Algols, and this was the 
aim of our study on TZ Eri. 
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Fig. 10. Mass-Luminosity relation for the components of the 
same Algol systems as those in Fig. 9. The zero-age main se- 
quence (ZAMS) and the Terminal- Age Main Sequence (TAMS) 
are drawn, according to Schaller et al (1992). Primaries are 
identified with filled symbols and secondaries with open sym- 
bols. TZ Eri is identified with plus symbols. Dots and squares 
refer to the systems having respectively under- and overlumi- 
nous primaries with respect to the ZAMS. 
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